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Abstract. Results of molecular dynamics simulations of fission reactions Na^Q — > 
Na^ + Na^ and Nalg — > 2iVaj^ are presented. Dependence of the fission barriers on 
isomer structure of the parent cluster is analyzed. It is demonstrated that the energy 
necessary for removing homothetic groups of atoms from the parent cluster is largely 
independent of the isomer form of the parent cluster. Importance of rearrangement 
of the cluster structure during the fission process is elucidated. This rearrangement 
may include transition to another isomer state of the parent cluster before actual 
separation of the daughter fragments begins and/or forming a "neck" between the 
separating fragments. 
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Fission of charged atomic clusters occurs when repulsive Coulomb forces, arising due 
to the excessive charge, overcome the electronic binding energy of the cluster 
This mechanism of the cluster fission is in a great deal similar to the nuclear fission 
phenomena. Experimentally, multiply charged metal clusters can be observed in the 
mass spectra when their size exceeds the critical size of stability, which depends on the 
metal species and cluster charge jUEHnj- 

We report the results of the ab initio molecular dynamics (MD) simulations of 
the fission processes NoIq Na^ + Na^ and Avals' ~* 2Na^ . Both symmetric and 
asymmetric fission channels are considered. We have investigated the parent cluster 
isomer dependence of the fission barrier for the reaction No^q — > Na\ + Na^ . To 
the best of our knowledge, a comparative study of fission barriers for various isomers 
by means of quantum chemistry methods has not been carried out before. Note that 
such a study is beyond the scope of simpler approaches which do not account for ionic 
structure of a cluster. 

We found that the direct separation barrier for the reaction Na\Q Na^ + Na^ 
has a weak dependence on the isomeric structure of the parent cluster. We note, however, 
that the groups of atoms to be removed from the parent cluster isomers must be chosen 
with care; one has to identify homothetic groups of atoms in each fissioning isomer. 
The weak dependence on the isomeric state of the parent Nal^ cluster implies that 
the particular ionic structure of the cluster is largely insignificant for the shape and 
height of the fission barrier. This is due to the fact that the maximum of fission barriers 
in considered cases are located at distances comparable or exceeding the sum of the 
resulting fragments radii. At such distances the interaction between the fragments, 
apart from the Coulombic repulsion, is mainly determined by the electronic properties 
rather than by the details of the ionic structure of the fragments. This is an important 
argument for justification of the jellium model approach to the description of the fission 
process of multiply charged metal clusters. 

We have demonstrated the importance of rearrangement of the cluster ionic 
structure during the fission process. The possibility of rearrangement of the cluster 
structure leads to the fact that direct fission of a cluster isomer in some cases may not 
be the energetically optimum path for the fission reaction. Alternatively, the reaction 
can go through transition to another isomer state of the parent cluster. This transition 
can occur in the first phase of the fission process, before separation of the fragments 
actually begins. We show that this is the case for the fission of and isomers of 
NoIq cluster. 

The rearrangement of ionic structure may be important also after the fragments 
began to separate. For Nafg — 2A^a9' reaction, two magic fragments Nag form a 
metastable transitional state in which the fragments are connected by a "neck". This 
"necking" allows for significant reduction in the height of the fission barrier. Note that 
the similar necking phenomenon is known for the nuclear fission process [7j. 

In our molecular dynamics simulations we utilize methods of density functional 



3 



theory (DFT). Within the DFT one has to solve the Kohn-Sham equations jS] 
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Ui + Vn + V^cj^i = eitPi, (1) 

where the first term corresponds to the kinetic energy of an electron from the Kohn- 
Sham reference system, Ui describes the attraction of the i^^ electron to the nuclei in 
the cluster, -01 is the electronic orbital, Vh is the Hartree part of the inter-electronic 
interaction, 

^^^^^ = / F^l^^' 

p(r) is the electron density, V^c is the local exchange-correlation potential defined as the 
functional derivative of the exchange-correlation energy functional 

Ko = (3) 

6p{r) 

where the exchange-correlation energy is partitioned into two parts, referred to as 
exchange and correlation parts: 

E^,[p] = E^{p) + E,{p). (4) 

Physically, these two terms correspond to same-spin and mixed-spin interactions, 
respectively. Both parts are functionals of the electron density, which can be of two 
distinct types: either a local functional depending on the electron density p only or a 
gradient-corrected functional depending on the electron density and its gradient Vp. 

There is a variety of exchange-correlation functionals in the literature. We have 
used the three-parameter Becke-type gradient-corrected exchange functional with the 
gradient-corrected correlation functional of Lee, Yang, and Parr (B3LYP) ^j. For the 
explicit form of this functional we refer to the original papers [TIH ITT| IT^ . The B3LYP 
functional has proved to be a reliable tool for studying the structure and properties 
of small metal clusters. It provides high accuracy at comparatively low computational 
costs. For a discussion and a comparison with other approaches, see [ini IHj . Note that 
the density of the parent cluster and two daughter fragments (including the overlapping 
region before scission point) almost does not change during the fission process (by 
analogy with the deformed jellium model, see ^3 for more details). This means 
that the B3LYP method works adequately for any fragment separation distances, d, 
during the fission process. 

The calculations have been carried out with the use of the GAUSSIAN 98 software 
package The 6-311G(d) and LANL2DZ basis sets of primitive Gaussian functions 
have been used to expand the cluster orbitals jTHj. The 6-311G(d) basis has been used 
for simulations involving Naf^ cluster. This basis set takes into account electrons from 
all atomic orbitals, so that the dynamics of all particles in the system is taken into 
account. For Nafg cluster we have used more numerically efficient LANL2DZ basis, for 
which valent atomic electrons move in an effective core potential (see details in [TH]). 

To simulate the fission process we start from the optimized geometry of a cluster 
(for details of the geometry optimization procedure see fH^ I14j) and choose the atoms 
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the resulting fragments would consist of. The atoms chosen for a smaller fragment 
are shifted from their locations in the parent cluster to a certain distance. Then, the 
multidimensional potential energy surface, its gradient and forces with respect to the 
molecular coordinates are calculated. These quantities specify the direction along the 
surface in which the energy decreases the most rapidly and provide information for the 
determination of the next step for the moving atoms. If the fragments are removed not 
far enough from each other then the attractive forces prevailed over the repulsive ones 
and the fragments stuck together forming the unified cluster again. In the opposite 
situation the repulsive forces dominate and the fragments drift away from each other. 
The dependence of the total energy of the system on the fragment separation distance 
forms the fission barrier. The aim of our simulations is to find the fission pathway 
corresponding to the minimum of the fission barrier. 

There are usually many stable isomers of a cluster with energies slightly exceeding 
the energy of the ground state isomer. In order to analyze the isomer dependence of the 
fission barrier in the reaction NoIq — > iVa^ + Na'^ we have picked two energetically low- 
lying isomers with the point symmetry groups 6*4^ and D^^ differing from the distorted 
Td point symmetry group of the ground state parent NoIq cluster. Three isomer states 
of the NuIq cluster are shown in figure 




Figure 1. Three isomers of Naf^ cluster. From left to right: the ground state isomer 
of distorted Td point symmetry group (total energy is -1622.7063 a.u.); an isomer of 
point symmetry group (total energy is -1622.6888 a.u., that exceeds the lowest 
energy state by 0.476 eV); an isomer of £'4^ point symmetry group (total energy is 
-1622.6860 a.u., that exceeds the lowest energy state by 0.553 eV). The homothetic 
group of three atoms marked by black color. 



In figure 121 we show fission barriers for separation three atoms from the C^^, ^4^, 
and Td isomers of the A^aio" cluster. In this figure zero level of energy is chosen for 
each parent isomer separately and corresponds to the minimum of total energy of that 
isomer. The initial distances between the centers of mass of two (future) fragments are 
finite so that the barriers do not start at the origin. 

The barriers for all three channels are close. The weak sensitivity of the fission 
barrier on the isomeric states of the reactants can be explained if one notices that the 
barrier maxima are located at distances comparable to or exceeding the sum of the 
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Figure 2. Fission barriers for separating the homothetic group of three atoms (marked 
by black cofor in figure^ from three isomers of No^q cluster derived from molecular 
dynamics simulations (direct Na^Q Na^ + Na'^ fission channel). The barriers 
plotted versus distance between the centers of mass of the fragments. Solid, dashed, 
and dashed-dotted lines correspond to distorted Td, C^y^ and D^d point symmetry 
groups isomers of the parent cluster, respectively. Energies are measured from the 
energy of the ground state of the corresponding isomers, i.e. we plot E — £'Td(C4i,. 134^)1 
where E is the total energy of the system and Ex^{Ci^..Did) ^'"'^ ground energies of 
the Td, Ciy and D4^d isomer states of the parent Naf^ cluster, respectively. 

resulting fragments radii, that is not far from the scission point. At such distances 
the interaction between the fragments, apart from Coulombic repulsion, is mainly 
determined by the electronic properties rather than by the details of the ionic structure 
of the fragments. This is an important argument for justification of the jellium model 
approach to the description of the fission process of multiply charged metal clusters. 

It is important to note that the barriers presented in figure El are calculated in 
assumption that fission occurs for the fixed (given) isomers. However, since C^v and 
D411 isomers are not the lowest energy states of NoIq system, there could be other 
processes competing with fission. One of such processes is rearrangement of the cluster 
structure. 

Rearrangement of the cluster structure during the fission process may significantly 
reduce the fission barrier. Such rearrangement may occur before the actual separation 
of the daughter fragments begins or after that. 

Fission of C^y and D^d isomers of NoIq cluster is an example of situation where 
rearrangement of the cluster structure takes place before the fragments start to separate. 
In figure El we show schematically the total energies of the NoIq T^, C^y and D^d isomers 
and barriers for the transitions between those states. It is seen from the figure that 
transition to the ground (Tii) state with subsequent fission into the Naf and Na^ 
fragments, A^a^g (C'4« or D^d) A^Qio (^rf) ~^ + ^^3"; (shown by solid lines) is 
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Figure 3. Energy levels of some states of the NoIq system and schematic barriers for 
transitions between these states. Energies are measured from the energy of the ground 
Td state of the Nalt cluster. 
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Figure 4. Fission barrier for Nalg 2Nd^ channel derived from molecular dynamics 
simulations as a function of distance between the centers of mass of the fragments. 
Energy is measured from the energy of the ground state of the Na\^ cluster. 
The arrow shows position of the meta-stable transitional state, see also figure This 
results is in a good agreement with the results of the jellium model |15[I16| . 



preferred path for fission of both C^^ and D^^ isomers of the Na^ cluster and requires 
only about 0.2 eV for C^y isomer and 0.26 eV for D^d isomer. In contrast, the direct 
fission process, A^a^Q (C4„ or D4d) — > Nay +Na3', (shown by dashed hues) requires about 
0.5 eV. We also show the barrier for the transition between the C4„ and D^^d isomers. 

Another example of cluster structure rearrangement in the fission process is the 
Nalg — > 2A^a9' reaction. The fission barrier for this reaction is shown in figure |3] If 



Figure 5. Rearrangement of the cluster structure during the fission process Nafl — > 
2Nag . From left to right: ground state of the parent cluster; "necking" between the 
two fragments leads to a meta stable intermediate state and significantly reduces the 
fission barrier height; two Nag fragments drifting away from each other. 



Table 1. Smirmary of the fission barrier heights (eV). 
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two fragments of the parent cluster were not allowed to adjust their ionic structure 
the fission barrier would be about 1 eV. Rearrangement of the cluster structure allows 
to reduce the fission barrier down to 0.31 eV. During the fission process the daughter 
fragments start to drift away from each other and a "neck" forms between the fragments. 
Formation of the "neck" results in a metastable transitional state. The geometry of this 
state, as well as the geometry of the parent cluster are shown in figure El 

In table ^ we have summarized our results for the fission barrier heights and 
compared them with the results of other molecular dynamics simulations and with the 
predictions of the jellium model. 

We have investigated two aspects of charged metal cluster fission process: 
dependence of the fission barrier on isomer state of the parent cluster and importance 
of rearrangement of the cluster ionic structure during the fission process. 

We found that for a consistent choice of the atoms removed from the cluster the 
fission barrier for the reaction Na\Q — > Na^ + Na^ has a weak dependence on the initial 
isomer structure of the parent cluster. This implies that the particular ionic structure 
of the cluster is largely insignificant for the height of the fission barrier. This is an 
important argument for justification of the jellium model approach to the description 
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of the fission process of multiply charged metal clusters. 

We have shown importance of rearrangement of the cluster ionic structure during 
the fission process. The fission reaction can go through transition to another isomer 
state of the parent cluster. This transition can occur before actual separation of the 
fragments begins and/or "neck" between the separating fragments is formed. In any 
case the resulting fission barrier can be significantly lower compared to the one for the 
direct fission path. 
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